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We carry out a model independent study of resonant plioton scattering ofT dark matter (DM) 
particles. The DM particle xi can feature an electric or magnetic transition dipole moment which 
couples it with photons and a heavier neutral particle X2- Resonant photon scattering then takes 
place at a special energy E^^^ set by the masses of xi and X2, with the width of the resonance 
set by the size of the transition dipole moment. We compute the constraints on the parameter 
space of the model from stellar energy losses, data from SN 1987A, the Lyman-a forest, Big Bang 
nucleosynthesis, electro-weak precision measurements and accelerator searches. We show that the 
velocity broadening of the resonance plays an essential role for the possibility of the detection of 
a spectral feature originating from resonant photon-DM scattering. Depending upon the particle 
setup and the DM surface mass density, the favored range of DM particle masses lies between tens 
of keV and a few MeV, while the resonant photon absorption energy is predicted to be between tens 
of keV and few GeV. 

PACS numbers: 98.80.Cq, 98.80.Es, 95.35+d, 12.60.Jv, 14.80.Ly 



I. INTRODUCTION 



One of the crucial properties of the dark matter (DM) 
is the feebleness of its coupling to the electromagnetic 
field. The early decoupling of DM from the baryon- 
photon fluid is a basic ingredient of the current picture 
of structure formation, and various direct DM detection 
experiments set stringent limits on the coupling of DM 
with ordinary matter. The phenomenological possibili- 
ties of a charged [l| , or a milli-charged 0, [3| DM species, 
or that of DM featuring an electric or magnetic dipole 
moment Q were considered in several recent studies, all 
pointing towards a severe suppression of any coupling of 
the DM with photons. Significant absorption or scat- 
tering of photons by DM appears to ruled out, perhaps 
implying that DM does not cast shadows. 

In this analysis we investigate the possibility that, 
while the typical scattering cross section of DM with 
photons is very small, photons with the right energy can 
resonantly scatter off DM particles. We show that this 
resonant scattering might result in peculiar absorption 
features in the spectrum of distant sources. This effect 
can occur if the extension of the standard model of parti- 
cle physics required to accommodate a (neutral) DM par- 
ticle candidate xi also encompasses (1) a second, heavier 
neutral particle X2 and (2) an electric and/or magnetic 
transition dipole moment which couples the electromag- 
netic field to xi and X2- We also assume, for definiteness, 
that xi and X2 are fermionic fields. In this setting, there 
exists a special photon energy El^^ where the scattering 
cross-section of photons by DM is resonantly enhanced 
to the unitarity limit. If the resonance is broad enough. 
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and the cross section and DM column number density 
are large enough, the spectrum of distant photon sources 
might in principle feature a series of absorption lines cor- 
responding to DM halos at different redshifts. If such 
anomalous absorption features exist, not only would they 
provide a smoking gun for the particle nature of DM, but 
they could also potentially give information about the 
distribution of DM in the Universe. 



II. THE MODEL 

We adopt here a completely model-independent set- 
ting, where we indicate with mi_2 the masses of Xi,2, 
and consider the effective interaction Lagrangian 
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In the rest frame of xi, the photon-DM scattering medi- 
ated through an s-channel X2 exchange [see Fig. [Ija)] is 
resonant at the photon energy 
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~ , the 7-DM scattering cross section can be 
approximated with the relativistic Breit-Wigner (BW) 
formula 
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where pcM indicates the modulus of the center-of-mass 
momentum, s is the center-of-mass energy squared, F^^ 
is the total decay width of X2, and F 



X2^Xl7 



is the de- 
cay width of X2 into x\ and a photon. The value of 
cr-vY, at E-y = El°^ saturates the unitarity limit provided 
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1. Under this assumption, even 



if xi s-iid X2 featured interactions different in their de- 
tailed microscopic nature from those described by Eq. ([T|) 
(such as a transition milli-charge, or fermion-sfermion 
loops in neutralino DM models), the maximal resonant 
7-DM scattering cross section would always be given by 
a-yxiiE^) for E-y « E!i^°^. From Eq. H]), we compute 
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In the remainder of this study, for conciseness, we shall 
denote 




X2 



(b) 



FIG. 1: (a): Feynman diagram associated to the resonant 
photon scattering cross section described by Eq. ^ ; the black 
square indicates the effective vertex of the lagrangian ((T)) . (b) : 
The xi — X2 loop diagram contributing to the vacuum polar- 
ization tensor (|15p. 



^5~) DM density. When the quantity r 1 in Eq. (fTU|) for 



and, in order to maximize the scattering rate of photons 
by DM, we will assume a model with -B^j-^-y « 1. All the 
quantities above can be trivially rephrased in terms of 
TO2 and of the two ratios R = mi/m2 and rj = TO2/M as 
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where F = F/m2. 

Let us now turn to the effects of the resonant scatter- 
ing of photons emitted by a distant source. The mean 
specific intensity at the observed frequency as seen by 
an observer at redshift zq from the direction ip is given 

by 



where 
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e is the emissivity per unit comoving volume, and Tcs is 
the effective opacity. The latter can be cast as 
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where p is the DM density. To get a numerical feeling 
of whether the resonant scattering of photons leads to a 
sizable effect, we define 



mi 
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where E^m indicates an effective DM surface density, as- 
sociated with the integral along the line of sight of the 



Ej = E^°^ we expect a significant absorption for photon 
energies E^ E^°'^. 

Once a photon from a background source scatters off an 
intervening DM particle, the flux from the source itself is 
attenuated as long as the photon is diffused into an angle 
larger than the angular resolution of the instrument. The 
kinematics of the process closely resembles that of the rel- 
ativistic Compton scattering [5] or Thompson scattering 
at lower energies. Roughly speaking, the relevant quan- 
tity can be cast as the fraction F of scattered photons 
which end up being scattered into an angle smaller than 
the instrumental angular resolution 0ar, over the total 
number of scattered photons. For an order of magnitude 
estimate it is easy to show that, apart from the details 
of the DM distribution geometry, F depends on the two 
variables Oar and x = E^/m^-^. Making simple assump- 
tions, we estimate the values of F for an instrument fea- 
turing an angular resolution of one degree, over the range 
10-2 < a; < 102, fall within lO^" < F < 5 x 10"^. We 
therefore can safely assume that if a photon scatters off 
DM, it is effectively lost {i.e. the flux of photons from the 
background source is effectively depleted by photon-DM 
scattering processes). 



III. CONSTRAINTS ON THE PARAMETER 
SPACE 

Since the lagrangian ^ effectively couples the elec- 
tromagnetic field to xi and X2, depending upon the size 
of the coupling and the mass of the two particles Xi,2, 
the constraints that apply to a milli-charged particle ip 
(e.g. neutrinos featuring a small electric charge e' ^ e 
[(^, or the paratons of Ref. 0) will also be relevant for 
the present setting. 

The parameter space of the model we consider here 
consists of the parameters m2, rj and R; for future conve- 
nience, we choose to represent the viable range of param- 
eters on the (rj, 7712) plane, at fixed representative values 
of i? = 0.01, 0.1 and 0.99 [Fig. [J. 

To translate the constraints from the milli-charged sce- 
nario in the present language, we need to compute the 
cross section cr(/^/^ — » X1X2), and compare it to the 
standard cr{f'^f^ — > ^ip) = e'^a{f^f~ e+e") cross 
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FIG. 2: Constraints from Lyman-Q, SN 1987 A, BBN, Elec- 
troweak precision and other collider data on the model under 
consideration here, on the (77, 7712) plane, at i? = 0.01, 0.1 and 
0.99. The parameter space regions ruled out lie below the 
curves in the figure. 



section, where e = 
f^(/^/" X1X2) 



e'/e. We find 
16a „ 2 



(11) 



l + (l-i?2)2^-2(l+i?2). 



A first simple astrophysical constraint on the model 
is based on avoiding excessive energy losses in stars that 
can produce X1X2 pairs by various reactions, in particular 
through plasma decay processes. The most stringent lim- 
its come from avoiding an unacceptable delay of helium 
ignition in low-mass red giants. The relevant energy scale 
in the process is the plasma frequency LOp k, 8.6 keV, and 
the limit applies, roughly, to masses TO2 < 2ujp/{\ + R), 
constraining 0] 

8 X 10"^^ < ^f?^ < 2 X 10"^^ [Red Giants]. (12) 

While the lower limit stems from the energy losses argu- 
ment, the upper limit comes from the requirement that 
the mean free path of xi is smaller than the physical size 
of the stellar core: if the xi particles get trapped, the 
impact on the stellar evolution through energy transfer 
would in any case be negligible compared to other mech- 
anisms 0,111 . 

At such low masses, however, constraints from large 
scale structure, and namely from Lyman-a forest data, 
on the smallest possible mass for the DM particle force 
R'm2 ^ 10 keV This bound corresponds to the left- 
most horizontal lines in Fig. O 

For a narrow range of effective 7 — Xi ~ X2 couplings, 
the cooling limit discussed above can be applied to the 



SN 1987A data. For a SN core plasma frequency ojp « 10 
MeV, values of TO2 ^ 2ujp/(l + R) can be ruled out in 
the range of couplings Q 

2 X 10^21 <Rrj'^<2x IQ-^'^ [SN 1987A]. (13) 

The limits from SN 1987A show up on the (77, 7712) plane 
as the rectangular regions of parameter space shown in 
Fig.H 

If X 1,2 reach thermal equilibrium in the Early Universe 
before big bang nucleosynthesis (BBN), they contribute 
to the energy density and thus to the expansion rate. 
Translating in the present language the constraints from 
BBN found in Ref. 0, if TO2 < 2 x (1 MeV)/(l + i?) then 
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The BBN limit corresponds to the central horizontal lines 
shown in Fig. [51 

As pointed out in , a strong constraint on the size of 
DM magnetic or electric dipole moments is related to the 
size of the photon transverse vacuum-polarization tensor 
[see Fig. a (b)]. 
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The strongest constraint derived from (jlSp comes from 
the effect of the running of the fine-structure constant, for 
momenta ranging up to the mass, on the relationship 
between mw, rnz and Gp. Using ([1]), we computed 
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The theoretically computed standard model values and 
the experimental inputs yield a limit on extra contribu- 
tions to the running of a, namely Ar < 0.003 at 95% 
C.L. [3. With particle masses 7771,7772 <C niz, Eq. pT]) 
reduces to 
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implying M > SAuiz for consistency with electroweak 
precision observables. The limits from Eq. (|17p rule out 
the region below the line labeled as "Vac. Pol. -EW preci- 
sion" in Fig. O 

Lastly, high energy accelerator experiments also put 
constraints on particles with an effective coupling to pho- 
tons. Such particles could have been seen in free quark 
searches Q, at the anomalous single photon (ASP) detec- 
tor at the SLAG storage ring PEP (designed to look 
for events in the form e+e~ ^7-1- weakly interacting 



4 




510 



FIG. 3: The photon-DM cross section at zero velocity (red 
dashed line) and with a velocity-dispersion cr^ = 820km/s 
(black line), for mi = 10 MeV, m2 = 100 MeV and F = 8.7 
eV. 



particles) and in beam-dump experiments from vector- 
meson decays and direct Drell-Yan production The 
combination of all accelerator constraints rules out the 
relatively massive and strongly coupled models lying be- 
low the upper-right curvy lines on the (ry, 7712) plane 
shown, for three values of R, in Fig. O 

This completes our discussion of the constraint on the 
parameter space of the model under consideration here: 
the viable parameter space, for a given R, lies above the 
lines shown in Fig. [51 while the portions of parameter 
space that are ruled out correspond to the regions of the 
plot below the various constraint lines. 



IV. RESONANT DM-PHOTON SCATTERING 

Since DM particles live in halos characterized by a ve- 
locity dispersion a^, which depending upon the mass of 
the DM halo can take values from roughly lOOkm/s to 
over 1000 km/s, the momentum distribution of the DM 
particles approximately follows a Maxwell-Boltzmann 
distribution. 
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An incoming photon will therefore scatter off DM par- 
ticles with the above momentum distribution, and the 
"ejffeciit'e" scattering cross section will be given by the 
following average: 
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where fi is the cosine of the incident DM-7 angle, and 
where the center of mass energy and momentum squared 
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FIG. 4: Curves of constant A2E/E!y°'' (i.e. the relative width 
of the DM-photon resonant scattering cross section) on the 
{av , R) plane. 
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The integral in Eq. (|2ip can be solved analytically, and 
we report the result in the Appendix. As a result of the 
averaging procedure of Eq. (j20p . the maximum of the 
effective cross section is no longer the peak value (T(i?"''), 
but will be a non-trivial combination of the latter, F^^ , 
cr„, and E!^°^. We illustrate an instance of the result of 
the broadening of the BW cross section in Fig. [3l 

Given an instrument with an energy resolution ^, de- 
fined as the relative energy resolution (i.e. the ratio of 
the energy resolution at the energy E over the energy E 
itself), we require the width A2E of the resonance (which 
we define, for convenience, to be the range of values of Ej 
where ct^^^ (E^y) > 10~^(t(£'^°*')) to be at least as large as 

To a good approximation, the solution to the equation 
A2E — ^xE^"^ is independent of ry, since cr^^^ (E^) oc F^ 



for E~, 



X2 

. Also, since A2E oc mi [see Eq. ^^], at 
fixed R and small F^^^ , the ratio A2E/El^'^ is independent 
of TO2 as well. We therefore plot, in Fig. [51 curves at 
constant values of A2E/E^°^ on the {av,R) plane. As 
clear from Eq. (fTOl) . the larger the value of the velocity 
dispersion a^, the larger A2E. From Eq. ([6]) we also 
understand that, as i? ^ 1, 0, explaining why 



arbitrarily large values of A2E/El^^ can be obtained for 
large R [see the upper part of Fig. [3]. 

How would the spectrum of a background source look 
like after photons have resonantly scattered off DM? We 
address this question in Fig. [51 We assume for definite- 
ness (our results don't critically depend upon the partic- 
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FIG. 5: The effect of photon-DM scattering on a putative 
power-law spectrum (see Eq. (|24|l ') for a source located behind 
(or at the center of) a cluster with features comparable to 
those of the Coma cluster. We assume = 10 MeV, 771^2 — 
100 MeV and four different values of the ^^resonance width" 
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ular spectral shape) a power-law spectrum of the form 

. (24) 

We consider a setup where m^^ — 10 MeV and m^^ — 
100 MeV, and as an example we focus on the case of a 
source located behind (or at the center of) a cluster with 
features similar to those of the Coma cluster. Making use 
of the estimates provided in Ref . [ll| , we consider a DM 
surface density (integrating the D05 DM profile ^12] along 
the direction of the center of the cluster, within one virial 
radius of the cluster center) of Sqm — 5 x 10^^ MeV/cm^. 
Also, we assume a velocity dispersion of (jy = 820 km/s. 
Notice that the redshift of the Coma cluster, z ~ 0.0231, 
is small enough that the effect of photon redshift on the 
shape and location of the absorption feature is completely 
negligible. Making use of these estimates, the effect on 
the background source spectrum depends entirely upon 
the value of T^.^ : for large values of the latter quantity 
the DM halo is opaque to photons with energies around 

^res 

We show in Fig.[5]how the spectrum defined in Eq. 
is affected by setups with various different values of T^^^ . 
For ^ 6 eV (orange line), the absorption is almost 
complete around i?"^. Smaller values of F^^ imply only 
a partial deformation of the spectrum, and a reduced en- 
ergy range where absorption effectively takes place. For 
< 1 eV the absorption feature would be almost in- 
visible. 

We summarize our results on the (77,7712) plane in 
Fig. [51 for the same reference values we employed in 
Fig. m i.e. i? = 0.1 and ay = 820 km/s. For this choice 
of parameters, /S.2E/E'^^^ k, 0.015. The area shaded 
in yellow at the bottom right of the plot is ruled out 
by the various constraints discussed in the previous sec- 



FIG. 6: The (7;, 7712) plane at i? = 0.1 and (t„ = 820 km/s. 
The region shaded in yellow is ruled out by the constraints 
discussed in Sec. IIIII The green dashed lines correspond to 
fixed values of the DM surface density required to have a sub- 
stantial absorption feature in the spectrum of distant sources, 
in units of S. = 10^" MeV/cm^ (blue solid line). 



tion. The green dashed lines correspond to fixed val- 
ues for Sdm such that r ^ 1 in Eq. (jlOp . in units of 
= lO^o MeV/cm^ (blue solid hne). For DM surface 
densities E^m ^ 0(10'^° MeV/cm^), absorption is possi- 
ble for DM particle masses toi = Rm2 « 0(0.3-50) MeV. 
The absorption feature, in this plot, is predicted accord- 
ing to Eq. ([2]) to occur at a photon energy El^'^ ~ 5?7i2; 
henceforth, in the range above, we predict El^^ « 0(15- 
2500) MeV. 

The analogue of Fig. [6] for other values of and R can 
be directly read out of our results shown in Fig. [¥] taking 
into account the constraints shown in Fig. [21 and the fact 
that the values of E^m such that t ~ 1 in Eq. (fTUj) scale 
approximately linearly with R. For instance, again for 
DM surface densities Sqm ^ 10"^° MeV/cm^, we would 
predict a DM particle mass 20 keV < mi < 50 MeV for 
R = 0.01, and around 1-50 MeV for R = 0.99. Analo- 
gously, the location of the absorption feature is predicted 
in the range S^'^'' « 1 MeV to 150 GeV for R = 0.01, and 
at £;:^°" w 10-500 keV for R = 0.99. In the present setup, 
therefore, for reasonable DM surface densities, the loca- 
tion of the absorption feature varies in a wide range of 
photon energies, from tens of keV up to several GeV. 



V. DISCUSSION 

Photons from background sources will in general pass 
through various DM halos at all intermediate redshifts. 
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resulting in a cumulative cosmological effect leading, in 
principle, to a broadening and modulation of the absorp- 
tion feature described above. In Ref. [13], for instance, an 
analogous computation was carried out for the monochro- 
matic photon emission from DM pair annihilations into 
two photons; the detailed setup here is however different 
as the effect depends linearly rather than quadratically 
on the dark matter density distribution. A similar cos- 
mological broadening was also discussed for the case of 
resonant vv — > high energy neutrino absorption e.g. 
in Ref. [31 ■ The spatial homogeneity of the cosmic neu- 
trino background results however in a completely differ- 
ent column density structure than in the present case. 
The detailed computation of this cumulative cosmologi- 
cal effect depends on several assumptions about the dis- 
tribution and nature of DM structures in the Universe, on 
the presence of DM clumps or other substructures [l5|, 
and on the assumed halo density profiles and velocity 
distributions We leave the detailed analysis of this 
effect to a future study. 

In passing we note that thermally-produced DM candi- 
dates with masses in the tens of keV to the MeV range, 
often referred to as warm DM candidates, exhibit po- 
tentially interesting features in structure formation, sup- 
pressing, through free-streaming, small scale structures, 
and partly alleviating the cusp problem of cold DM mod- 
els (see e.g. Ref. ^] and references therein). Depending 
on the details of the particle physics model constraints 
on such warm DM candidates might be used to constrain 
our scenario. 

Closing the photon line in Fig. [1] (a) into a loop gen- 
erates radiative corrections to mi. If the latter are too 
large, the values we employed must be corrected accord- 
ingly, and small values of mi might not be theoretically 
allowed. We can estimate the size of these corrections as 



5mi 1 m\ R^rf 
I67r2 TO2 M2 ~ 167r2 



(25) 



Radiative corrections are therefore smaller than 10^^ 
provided 77 < 1, a condition which is always widely sat- 
isfied in the parameter space under consideration here. 
A mass mixing term would also be generated by the in- 
teraction responsible for the effective lagrangian ([T]); in 
principle, one should then rotate Eq. ^ to the proper 
mass-eigenstate basis. However, the relative size of the 
induced xi ~ X2 mixing is also very suppressed, as it 
roughly scales as Rrf, and can be thus safely neglected 
here. 

In the scenario we are discussing here, the xi particles 
can also pair annihilate into two photons through a X2 t- 
or u- channel exchange. The resulting cross section can 
be estimated as 
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Pair annihilation of xi's into photons can a priori be 
the process through which DM annihilates in the early 



Universe and potentially this could thermally produce 
the amount of DM inferred in the current cosmological 
standard model. In the range of couplings and masses we 
obtain here, the above mentioned annihilation channel is 
insufficient to produce a large enough pair annihilation 
rate in the Early Universe in order to get the required 
DM abundance 17^^ sa f^cDM- Other channels, otherwise 
irrelevant for the present discussion, and compatible with 
the present setting, can however contribute to give the xi 
particles the right pair annihilation rate 

The same diagram discussed above, and the same pair 
annihilation cross section, intervene in the pair annihi- 
lation rate of xi's today into monochromatic photons of 
energy E-y^ = to^j^ — Rm2. The flux of photons per 
unit solid angle from monochromatic pair annihilations 
of xi's can be written as 



R^j]"^ j{e,An) 

rrvk Att 
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where, in this instance, the quantity J refers to the fol- 
lowing line-of-sight integral along the direction 9 aver- 
aged over the solid angle Afi 
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An 



AO 



dl p^{l). 



(28) 



When mi < m2 ^ rnz, we can derive an upper limit to 
the monochromatic photon flux which is independent of 
mi ^2, namely 



A17) < 0.3i?2 



J(6>,A0) 
GeV cm-= 



cm"2 s^^sr-i (29) 



Taking an angular region AO, ~ sr in the direc- 

tion of the galactic center the range of values which J 
can take for various viable DM halo models is 10"^- 
10^2 GcV* cm^"^. This means that one expects a 
flux of monochromatic photons in the range (I0~^- 
10^'^)i?2 cm^2 s"-'^sr^^. The diffuse gamma-ray flux in 
the Galactic center region as measured by COMPTEL 
and EGRET ^ is at the level of 0.01 cm-2 s-^sr"! at a 
gamma-ray energy of 1-3 MeV. Extrapolating to smaller 
energies we expect an even larger flux at energies around 
or smaller than 100 keV. This makes it extremely hard to 
reconstruct a would be annihilation signal from the galac- 
tic background. Dedicated searches for line emissions 
show that instruments such as INTEGRAL-SPI also fail 
to achieve the sensitivity required here . On the other 
hand, this also means that the class of models discussed 
above is not currently constrained by monochromatic 
photon emissions. Furthermore, observations of objects 
where the diffuse gamma-ray background is expected to 
be suppressed, such as the nearby dwarf galaxies [l9| , can 
potentially lead to constraints or even to the detection of 
the monochromatic emission line predicted here. 

If, as we describe here, photons scatter off DM at sig- 
nificant rates, one might also expect other associated fea- 
tures besides the absorption lines and the monochromatic 
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emissions described above. Scattering off DM might 
generate an effective "index of refractiori^ in the pho- 
ton propagation, possibly inducing e.g. time delays in 
transient sources at different frequencies, or frequency- 
dependent distortions of the photon paths for steady 
sources. A detailed discussion of these effects lies, how- 
ever, beyond the scopes of the present analysis. 

Neutralino DM in the context of the minimal super- 
symmetric extension of the Standard Model (MSSM) can 
in principle produce an effective lagrangian setup as that 
in Eq. ([1]), for instance through fermion-sfermion loops 
coupling two different neutralinos xi and Xi, i > 1. From 
the discussion above, however, it is clear that supersym- 
metric DM cannot produce any sizable photon absorp- 
tion. First, the lightest supersymmetric particle (LSP) 
in any viable low energy supersymmetry setup is typi- 
cally heavier than at least a few GeV (for exceptions, 
e.g. in the next-to-MSSM, see [23|)- This implies, as can 
be read off Fig. [HI very large values of E^m to get r 1 
in Eq. (fTO|) . Secondly, the assumptions we made at the 
beginning that -B^it ~ 1 does not hold in general in the 
MSSM: the radiative X2 Xi7 decay can be the domi- 
nant mode only in restricted regions of parameter space, 
e.g. when phase space suppresses other three-body de- 
cay modes. The resulting effective 77, in the notation set 
above, is in any case limited from above by 



167r2 



< 10" 



Requiring F-^ 



X2 



implies m 



X2 



(30) 



and typi- 



cally R > 0.995. Furthermore, since in the MSSM when 
two neutralinos are quasi degenerate the lightest chargino 
is also quasi degenerate with them, LEP2 limits on the 
chargino mass Q force m^^ > 100 GeV. R > 0.995 also 
implies Tj^^ ^ 10"^" GeV. These values for the model 
parameters imply (a) small relative widths and (6) too 
large DM surface densities for the absorption feature to 
be detectable. Relaxing the requirement that B^-^^ « 1 
would not help anyway, since the cross section ^ re- 
ceives the large suppression factor T^2->xi'y/^X2J the 
photon absorption process is again suppressed. 

One can envision, however, various particle physics 
scenarios where the phenomenology described above can 
take place. For instance, a concrete particle physics setup 
which can explain at once the DM abundance, neutrino 
masses and mixing, the baryon asymmetry of the Uni- 
verse and, potentially, inflation, is the so-called z^MSM 
[2l| . or one of its extensions These models feature 
a light quasi-stable sterile neutrino with a mass in the 
tens of keV (21)] up to 0(10) MeV [22| range, and heavier 
Majorana neutrinos with a mass at the GeV scale. Ex- 
tending this class of models with an effective interaction 
of the form of our Eq. ([T]) gives rise to the phenomenology 
described above and, hence, to possible resonant photon 
scattering. 



VI. CONCLUSIONS 

We have shown that photons can, in principle, reso- 
nantly scatter off DM, through an effective lagrangian 
featuring a dipole transition moment coupling photons, 
the DM particle xi and a heavier neutral particle X2- We 
discussed the constraints on the model from stellar en- 
ergy losses, data from SN 1987A, the Lyman-a forest, 
Big Bang nucleosynthesis, electro-weak precision mea- 
surements and accelerator searches. The effective reso- 
nant "absorption" cross section is broadened by the effect 
of the momentum distribution of DM particles in DM ha- 
los. We showed that DM particles in the tens of keV to 
a few MeV range can lead to resonant photon scatter- 
ing (resulting in absorption lines which can lie between 
tens of keV up to tens of GeV) provided the DM surface 
mass density is at least of 0(10^® MeV/cm^). We also 
pointed out that typical supersymmetric DM (the weak- 
scale neutralino) does not cast any shadows (i.e. it does 
not "absorb" photons), while photon absorption can take 
place in other particle physics setups which can explain 
various pieces of physics beyond the standard model. 
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APPENDIX: THE ANGULAR INTEGRAL (a)^ 

The angular integral (a)^ can be computed analyti- 
cally with the result 



Mm,T,,r X (A.l) 



£2 [(Am2)2+r7i2r2j 
|2 + ^[c^(A-/2) + c,g]| 



where 



Am2 = to2 - to2 



c/ = 1 



2mfAm'^ 



(Am2)2 + mlTl^^ 



(A.2) 



(A.3) 
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2Am2 



X2 



(Am2) 



(A.4) 



and 



{2ml ~ Am2)m2r^ 



(Am2) 



2^2 



m2r2^ 



/i=ln 



(2£;^v/p^+m2-2i;^p)2 



/2=ln 



(2i?^V^ 



2^;^,^ - Ato2)2 



(A.5) 



m2r2^j 



(A.6) 



arctan 



arctan 



Am2 + 2E, 



(p- \fn^ 



Am2 - 2E^ + + p2 



(A.7) 
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